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Abstract 



This paper is an extension of Fujii et al ( quant-ph/0307066 ) and in this one we again 



treat a model of atom with n energy levels interacting with n(n-l)/2 external làser fields, 
which is a natural extension of usual two level system. Then the rotating wave approxi- 
mation (RWA) is assumed from the beginning. 

To solve the Schròdinger equation we set the consistency condition in our terminol- 
ogy and reduce it to a matrix equation with symmetric matrix Q consisting of coupling 
constants. 

However, to calculate exp(— itQ) explicitly is not easy In the case of three and four 
level systems we determine it in a complete manner, so our model in these levels becomes 
realistic. 

In last, we make a comment on Cavity QED quantum computation based on three 
energy levels of àtoms as a forthcoming target. 



*E-mail address : fujii@yokohama-cu.ac.jp 
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1 Introduction 



The purpose of this paper is to develop a useful method applicable to Quantum Computation. 
As an easy introduction to it see for example 0, [Sj- 

Quantum Computation is in a usual understanding based on qubits which are based on two 
level system (two energy levels or fundamental spins of àtoms), See jl], [E], [Hj as for general 
theory of two level system. 

In a realistic image of Quantum Computer we need at least one hundred àtoms. However, 
we meet a very severe problem called Decoherence which may destroy a superposition of 
quantum states in the process of unitary evolution of the system. See for example [Zj or recent 
jHj as an introduction. At the present time it is not easy to control Decoherence. 

By the way, an atom has in general infinitely many energy levels, while in a qubit method 
only two energy levels are used, see for example 8J as an introduction of two level approximation. 
We should use this possibility to reduce a number of àtoms. Since it is not realistic to take all 
energy levels into consideration at the same time we use n energy levels from the ground state, 
which is in general called qudit theory. See for example [TT], [T2], jHj and [THj . 

In the paper ^Hj we considered a model of atom with n energy levels interacting with n(n- 
l)/2 external làser fields, which is a natural extension of usual two level system. The rotating 
wave approximation (RWA) is assumed from the beginning. 

In the model it is assumed that all coupling constants regarding interactions of atom with 
different làser fields are equal, which is nothing but an approximation theory. To construct 
a realistic model we again treat a "full" model of atom with n energy levels interacting with 
n(n-l)/2 làser fields. 

To solve the Schròdinger equation we set the consistency condition in our terminology 
and reduce it to a matrix equation with symmetric real matrix Q consisting of all coupling 
constants. However, it is not easy to calculate exp(— üQ) explicitly. 

Therefore we restrict to special cases. In the case of three and four level systems we deter- 
mine it in a complete manner, so they become realistic. See also |17j . 

In last, we make a comment on Cavity QED quantum computation based on three energy 
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levels of àtoms that is our forthcoming target. 



2 Review on Two Level System 



Let us review a two level system within our necessity and its Rabi oscillation. Let {ai, 02, 03} 
be famous Pauli matrices 
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Let us consider an atom with two energy levels E Q and E\ (Ei > E ) as two level approxi- 
mation. Its Hamiltonian is in the diagonal form given as 

/ 
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This is rewritten as 



Hq = E 
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where A = Ex — E is the energy difference. Since we usually take no interest in constant terms, 
we can set 

#o = y (la - 0-3) . (3) 

We consider an atom with two energy levels which interacts with external (periòdic) field 
with g cos(ujt + <p). In the following we set h — 1 for simplicity. The Hamiltonian in the dipole 
approximation is given by 



H = H + gcos(iüt + 4>)ai = y (1 2 - 03) + gcos(ut + (j))ai, 



(4) 



where u is the frequency of the external field, g the coupling constant between the external 
field and the atom. This model is complicated enough to solve, see jH], [Tü] . [TE], [Tïïj . 



In the following we set = for simplicity and assume the rotating wave approximation 
(which neglects the fast oscillating terms), namely 



cos(a;í) = ~(e itJÍ + e~ iuJt ) = -e iuJt (l + e~ 2iuJt ) « -e iu}t 



and 
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therefore the Hamiltonian is given by 
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by the redefinition of g (g/2 — ► g). It is explicitly 
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We would like to solve the Schròdinger equation 



d 



i—q = Hm. 

dt 



For that purpose let us decompose H in (jïïj) into 
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so if we set 
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then it is not difhcult to see 
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which is easily solved. For simplicity we set the resonance condition 

A = tu, 



(6) 
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then the solution of (fTÜj) is 



<]>(/) =exp{-igt\ ° 1 j }<I>(0) 
1 



/ 



cos(gt) — i sin(gí) 



i — isin(gt) cos(gt) 
As a result, the solution of the equation (J7J) is given as 



$(0). 



*(f) 
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d>(0) (12) 



by (jÜJ). If we choose $(0) = (1, 0) T as an initial condition, then 



cos(gt) 



_2 e lujt sm(gt) 

This is a well-known model of the Rabi oscillation (or coherent oscillation). 



(13) 



3 General Theory of N Level System 

In general, an atom has an infmitely many energy levels, however it is not realistic to consider 
all of them at the same time. Therefore we take only n energy levels from the ground state 
into consideration (E n _i > ■ ■ • > E\ > E ). Then from the lesson in the preceding section the 
energy Hamiltonian can be written as 







Ai 



(14) 



where Aj = Ej — E for 1 < j < n — 1. For this atom we consider the interaction with n(n — 1)/2 
independent external làser fields corresponding to every energy difference (Ej — Ei for j > i). 
For example see the following figure for the case of n = 4 : 
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Fig.l Atom with four energy levels and general action by làser fields 



Then the interaction term assuming the RWA from the beginning is given as 
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where {#0/3} are coupling constants, so the Hamiltonian is 

H = H + V. 

Here we set 

= ^01) ^2 = ^12) • • • j ^n-2 = ^n-3,n-2) w n-l = w n-2,n-l 

for simplicity. We would like to solve the Schròdinger equation 

dt 
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Similarly in (JHJ) let us decompose H : For 
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it is easy to see 
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where 

for j — i > 2. 
By setting 

it is not difncult to see 

where 



— k>ij — (c^i+l + ^i+2 + ' ' ' + ^jy 



(20) 
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At this stage we take the resonance conditions 
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(22) 



to make the situation simpler. 

Moreover, we consider a very special case : namely, in (}2T]) we set 



tij = 



&ij = ^>i+i + ^i+2 + • • • + tjjj = Ej — Ei 



(23) 



for all j—i>2 (see for example the figure 1 once more). We call this a consistency condition. 
Then Hu becomes a constant symmetric matrix ! 
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It is easy to solve the equation (J20J) with constant Q 
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whose solution is formally given by 

f (í) = exp(-itQ)* (0). 
As a result we have a general solution 

*(í) = í/exp(-iíQ)#(0). 



(26) 



(27) 



Therefore the problem left is to calculate exp(— itQ) explicitly, which is however a very 
hard task. In the following, let us treat the special cases n = 3 and n = 4 because the formula 
in becomes helpful in the cases. 



4 Exact Solution in Three Level System 

In this section we consider the case of n = 3 and calculate exp(— itQ) in a complete manner. 

For simplicity we set u m = uj 1 , uj 02 = uj 3 , u 12 = uj 2 and g 01 = g 1 , g 02 = flfc, 9ví = 92- See the 
following figure. 
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Fig.2 Atom with three energy levels and general action by làser fields 
Now the matrix Q in (J21|) is in this case 
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(28) 



0= |A1 3 -Q| 



A 3 - (gï + g\ + gfj A - 2 gi g 2 g 3 



so we look for eigenvalues of Q. 
From 

A -0i -g 3 
-gi A -g 2 
-9a ~92 A 

the Cardano formula (see for example [21]) gives three real solutions 

Ai = a + + a_, A2 = u 2 a + + aa^, A3 = aa + + cr 2 o;_ 
where a = e 27 ™/ 3 and 



<X± = { 9v9-29:\ +- \UÍ\9m 
From [20] we have finally 
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(A 2 - Ax)(A 3 - Ai) (Ai - A 2 )(A 3 - A 2 ) (Ai - A 3 )(A 2 - A 3 ) 



-itXi 



-ÜX2 



-it\z 



+ 



(A 2 -Ai)(A 3 -Ai) (Ai-A 2 )(A 3 -A 2 ) (Ai - A 3 )(A 2 - A 3 ) 



(29) 



(30) 

(31) 
(32) 
(33) 



5 Exact Solution in Four Level System 

In this section we consider the case of n = 4 and calculate exp(— itQ) in a complete manner. 

For simplicity we also set u m = wi, ^02 = ^4, ^03 = <^6, ^12 = ^2, ^13 = &5, ^23 = ^3 and 
9oi = 9i, 902 = 94, 903 = 9e, 9u = 92, 9\3 = 9ò, 923 = 93- See the figure 1 once more. 

Now the matrix Q in (J21j) is i n this case 

^0 gi g 4 g e ^ 



Q 



9i g 2 g 5 
94 g2 g 3 
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(34) 



so we look for the eigenvalues of Q : 



= |A1 4 -Q| 



A -gi -94 -g e 
-g\ A -5-2 -05 
~9a ~92 A -5-3 
-0e -gs -93 A 

= A 4 - (0? + g\ + g\ + g\ + g\ + #j) A 2 - 2 (010204 + 010506 + 020305 + 03040è) A 

+ g\g\ + g\gl + g\g\ - 20102030è - 201030405 - 20204050è 

= A 4 + pA 2 + q\ + r. 



The Euler formula (see |21j ) gives four real solutions 

Ai = a + (3 + 7, A 2 = ara + a 2 /3 + a-y, A 3 = a 2 a + (3 + o 2 ^, A 4 = era + ex 2 /? + cr 3 7 (35) 

with very complicated terms a, (3 and 7 given below (corresponding to a+, a_ in the preceding 
section) and a = e 27U//4 = i. For the algebraic equation with three degrees 

64S 3 + 32pB 2 - 4(4r - p 2 )5 - q 2 = 0, 



we set (3 = y B for the largest real solution B (this can be obtained by using the Cardano 
formula again). a and 7 are given as primitive solutions of equations 
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From [2Ü] we have finally 



exp(-iíQ) = / (í)l 4 + /i(í)Q + f 2 (t)Q 2 + Ht)Q 3 
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6 Approximate Solution of N Level System 

Since we are interested in the whole of exp(— itQ) in (}4*Sj) we take an approximation to the 
equation. Namely, we assume that all coupling constants are equal : g = g a p for < a, (3 < 
n — 1. Then 
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1 1 ■ • 1 1 

1 1 • • 1 1 
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/ 



1 1 1 • • 1 
1 1 1 • • 1 
and it is not difficult to calculate exp(-itgR) explicitly, [Tïï] . 
If we define 

|1) = (1,1,···,1,1) T , 
then it is easy to see R = |1)(1| — l n , so 



exp(-ügR) = e i9í exp(-^í|l)(l|). 
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Since (1|1) = n, 



so that 



l)(l|) fc = ^|l)(l|, 
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Therefore we obtain 

exp(-üQ) = exp(-ügR) = e 



exp(-ingt) — 1 



n 



- / ln + exp(-m^)-l 
n 



(42) 



(43) 



7 Discussion 

In this paper we developed a general theory of the n level system of atom by assuming the 
rotating wave approximation and reduce the Schròdinger equation to the simple matrix equation 
with matrix consisting of coupling constants under the consistency condition. 
Moreover, we solved the equation in a perfect manner in the three and four level systems. 

By the way, to assume the rotating wave approximation is a bit weak in the theory, so that 
we must study the general equation for example in the case of three level one like 

( 0i C0S(CJ1Í + 0i) #3 C0S(CJ 3 Í + 3 ) ^ ( lT ' 

d 



'dt 



^2 



#1 



(44) 



#i cos(cüií + 0i) g 3 cos(uj 3 t + 3 ) 

gi cos(cüit + 0i) Ai ^2 cos(cj 2 í + 02) 

^ 5f 3 cos(cj 3 t + 3 ) 02 cos(cj 2 í + 02) A 2 J 

However, it is almost impossible to solve this at the present time. Let us leave it for young 
researchers in Quantum Òptics or Mathematical Physics as a challenging task. 

Here we state our motivation once more. We would like to construct a realistic model of 
quantum computation with three level system. Its candidate is a quantum computation based 
on Cavity QED making use of three energy levels of àtoms. See the following figure. 
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Fig.3 A general setting for a quantum computation based on Cavity QED 



Let us add a short explanation to this figure. The dotted line means a single photon 
inserted in the cavity and all curves mean external làser fields (which are treated as classical 
ones) subjected to ultracold-atoms trapped linearly in it. See in detail [23 and [21] in which 
the Cavity QED quantum computation with two level system was "completed" . 
This is our forthcoming target ! 
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Appendix Diagonalization of Q in the Three Level 

In the case of three level system let us calculate exp(— itQ) by use of the diagonalization 
method which is more popular. 

In ()29|) we had the three eigenvalues {Ai,À2,A3}. Therefore the normalized eigenvectors 
{|<\j) | 1 < 3 ' < 3} corresponding to these are given as 




(45) 
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This derivation is given in the latter half. Then the matrix defined by 



0=(|A 1 ),|A 2 ),|A 3 )) 



makes Q diagonal like 
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T . 



Therefore we obtain the desired form 
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It is assumed g^ = g\ in the right hand side. 



Let us calculate a special case #3 = 0. From 
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9i 



\ 



V 2( ~9Ï+9 2 2 ) 

1 
V2 

92 



o 



10} 



92 



\ 







-91 



gl + gl) 



91 



V 2 Ï9Ï+9l) 
1 

V2 
92 



#i 2 + #l),|0),|^ + ^ 2 ) 



a straightforward calculation leads us to 



/ g 2 COS(ty/g 2 +g 2 )+g 2 . giSÍn(f y/gj+g 2 ) gig2COS(ty/ g 2 +g 2 )- gi g 2 \ 



exp(-itQ) 



g'í+g'í 



Vaï+92 



9Í+9Ï 



_. gl sin(^Fg) cog(í r^—^ _. gggm(Vg^|) 



y/sj+gi 



gig2COS(ty/gl+gZ)- gi g2 . g 2 Ún{t^/ gl+g 2 2 ) g 2 2 COS{t^/g 2 +g 2 )+g 2 

9Ï+9Ï ^ v 7 ^^! 



1 In the calculation we must choose signs of the complex roots in 1)45(1 carefully 
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See §3 in [TT] . 



Derivation of (1451) 



For Q in (J2%|) we consider the equations 



Qxj = ÀjXj for 1 < j ' < 3. 



Namely, 
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9iV + 93^ = XjX 
gix + g 2 z = Xjy 
93X + g 2 y = Xjz 



We solve the above equations with respect to z. The result is 
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where we have used the characteristic equation 
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Next let us determine the normalization. 
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Now we want to remove the term Ag^g^Xj. By using the characteristic equation 2g 1 g 2 fi , 3 
Af - (9Ï + 92+ 9Ï) X J, we have 

The right hand side of © = 3Aj - 3<7 2 A 2 - (gf + g\ + gl)\) + g\{g\ + £ 2 + g\) 
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from (jïïTJ). At this stage it is easy to conjecture 



Xj93 + 9i92 



\]~9l 



A 2 - gl 



Xj92 + 9x9^ 



2 - n 2 
j 9l 



A| - ^ 



In fact, 



(55) 



Aj ^3 + 9i92 

aF^! 



(Aj#3 + 



\ 



# 2 Af + 19\929?,\ + 9Í9Í 



where we have used the equation 



\ 2 n 2 

A i - 9i 



\ 



A 2 - gl 



\ 2 n 2 

A 7 - — g2i 



g\X) + ïgigm ^i + sUsü = & ó + Aj - {g\ + g\ + g\)\) + g\g\ 

= ^-{gl + gl)^ + glgl 
= tf-g 2 m-gl)- 



Similarly, we obtain 



Aj ff2 + gig-í 



A? 



gl 
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From (|55|) we finally obtain 



1 



- g\ - gl 



à 



V 



'\ 2 


-ai 


% 2 


-à 



^ 2 y 



v 





A 2 - 

3 






ni) 

3A 2 


-sí 




~9Í 




A 2 - 
j 


-4 






-sí 


-9Ï 


-9Í 




A 2 - 
j 






3A Í 


-9Ï 


-al 


SÍ 



A comment is in order. y/t means the complex square root defined by y/t = e , argí / 2 
in our case we must choose its sign (for example, a/4 = 2 or — 2) carefully. 



so 
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